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O ' Abstract 

. The hadronization of an expanding partonic fireball is studied within the Parton-Hadron-Strings 

I Dynamics (PHSD) approach which is based on a dynamical quasiparticle model (DQPM) matched 

" to reproduce lattice QCD results in thermodynamic equilibrium. Apart from strong parton inter- 

. actions the expansion and development of collective flow is found to be driven by strong gradients 

I in the parton mean-fields. An analysis of the elliptic flow V2 demonstrates a linear correlation with 

^ ' the spatial eccentricity e as in case of ideal hydrodynamics. The hadronization occurs by quark- 

^H^' antiquark fusion or 3 quark/3 antiquark recombination which is described by covariant transition 

! rates. Since the dynamical quarks become very massive, the formed resonant 'pre-hadronic' color- 

^ I dipole states (gg or ggg) are of high invariant mass, too, and sequentially decay to the ground- 

■ state meson and baryon octets increasing the total entropy. This solves the entropy problem in 

, hadronization in a natural way. Hadronic particle ratios turn out to be in line with those from a 

I grandcanonical partition function at temperature T ~ 170 MeV. 

d ; 

oo 1 Introduction 

o : 

>- ', The 'Big Bang' scenario implies that in the first micro-seconds of the universe the entire state has 
^ ', emerged from a partonic system of quarks, antiquarks and gluons - a quark-gluon plasma (QGP) - to 
color neutral hadronic matter consisting of interacting hadronic states (and resonances) in which the 
■ - - partonic degrees of freedom are confined. The nature of confinement and the dynamics of this phase 
transition has motivated a large community for several decades and is still an outstanding question of 
todays physics. Early concepts of the QGP were guided by the idea of a weakly interacting system of 
partons which might be described by perturbative QCD (pQCD). However, experimental observations 
at the Relativistic-Heavy-Ion Collider (RHIC) indicated that the new medium created in ultrarelativistic 
Au+Au collisions was interacting more strongly than hadronic matter and consequently this concept 
had to be given up. Moreover, in line with theoretical studies in Refs. [1, 2, 3] the medium showed 
phenomena of an almost perfect liquid of partons [4] as extracted from the strong radial expansion and 
elliptic flow of hadrons [4] . 

The question about the properties of this (nonperturbative) QGP liquid is discussed controversially 
in the literature and dynamical concepts describing the formation of color neutral hadrons from partons 
are scarce [5, 6, 7, 8]. A fundamental issue for hadronization models is the conservation of 4-momentum 
as well as the entropy problem because by fusion/coalescence of massless (or low constituent mass) 
partons to color neutral bound states of low invariant mass (e.g. pions) the number of degrees of 
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freedom and thus the total entropy is reduced in the hadronization process [6, 7]. This problem - a 
violation of the second law of thermodynamics as well as of the conservation of four-momentum and 
flavor currents - definitely needs a sound dynamical solution. 

A consistent dynamical approach - valid also for strongly interacting systems - can be formulated 
on the basis of Kadanoff-Baym (KB) equations [9, 10] or off-shell transport equations in phase-space 
representation, respectively [10, 11]. In the KB theory the field quanta arc described in terms of 
propagators with complex selfenergies. Whereas the real part of the selfenergies can be related to 
mean-field potentials, the imaginary parts provide information about the lifetime and/or reaction rates 
of time-like 'particles' [3] . Once the proper (complex) selfenergies of the degrees of freedom are known 
the time evolution of the system is fully governed by off-shell transport equations (as described in Refs. 
[10, 11]). 

The determination/extraction of complex selfenergies for the partonic degrees of freedom has been 
performed in Refs. [3, 12, 13] by fitting lattice QCD (IQCD) 'data' within the Dynamical QuasiParticle 
Model (DQPM). In fact, the DQPM allows for a simple and transparent interpretation of lattice QCD 
results for thermodynamic quantities as well as correlators and leads to effective strongly interacting 
partonic quasiparticles with broad spectral functions. 



2 The PHSD approach 

The Parton-Hadron-String-Dynamics (PHSD) approach is a microscopic covariant transport model 
that incorporates effective partonic as well as hadronic degrees of freedom and involves a dynamical 
description of the hadronization process from partonic to hadronic matter [14]. Whereas the hadronic 
part is essentially equivalent to the conventional Hadron-Strings-Dynamics (HSD) approach [15] the 
partonic dynamics is based on the Dynamical QuasiParticle Model (DQPM) [12, 13] which describes 
QCD properties in terms of single-particle Green's functions (in the sense of a two-particle irreducible 
(2PI) approach). 

We briefiy recall the basic assumptions of the DQPM: Following Ref. [16] the dynamical quasiparticle 
mass (for gluons and quarks) is assumed to be given by the thermal mass in the asymptotic high- 
momentum regime, which is proportional to the coupling g{T/Tc) and the temperature T with a running 
coupling (squared), 

g (T/n) = ^11^^ _ 2Nf) ln[A2(T/T, - TjT.y^ ' 

Here Nc = 3 stands for the number of colors while Nf denotes the number of flavors. The parameters 
controlling the infrared enhancement of the coupling A = 2.42 and Tg = 0.46Tc have been fitted in [16] 
to IQCD results for the entropy density s{T). An almost perfect reproduction of the energy density 
e{T) and the pressure -P(T) from IQCD is achieved as well. 

The width for gluons and quarks (for vanishing chemical potential Hg) is adopted in the form 



^g^T , f2c\ g'^T , /2c\ 



where c— 14.4 (from Ref. [3]) is related to a magnetic cut-off. We stress that a non- vanishing width 7 is 
the central difference of the DQPM to conventional quasiparticle models [17]. It influence is essentially 
seen in correlation functions as e.g. in the stationary limit of the correlation function in the off-diagonal 
elements of the energy-momentum tensor T^'' which defines the shear viscosity rj of the medium [3]. 
Here a sizable width is mandatory to obtain a small ratio in the shear viscosity to entropy density rj/s. 
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In line with [16] the parton spectral functions thus are no longer S— functions in the invariant mass 
squared but taken as 




(3) 



separately for quarks and gluons (j = q^q^g)- With the convention E'^{p) = + Mj — 7^, the 
parameters M| and 7^ are directly related to the real and imaginary parts of the retarded self-energy, 
e.g. Uj = Mj — 2i'jjU. With the spectral functions fixed by Eqs. (l)-(3) the total energy density in the 
DQPM (at vanishing quark chemical potential) can be evaluated as 

where ub and np denote the Bose and Fermi functions, respectively. The number of transverse gluonic 
degrees of freedom is = 16 while the fermic degrees of freedom amount to dg — ANcNf — 36 in 
case of three flavors (Nf—S). The pressure P then may be obtained by integrating the differential 
thermodynamic relation 

P-T^- _rOO (5) 

with the entropy density s given by 

dP TOO + p 

This approach is thermodynamically consistent and represents a 2PI approximation to hot QCD (once 
the free parameters in (1) and (2) are fitted to lattice QCD results as in Refs. [3, 12, 13]). 

As outlined in detail in Refs. [12, 13] the energy density functional (4) can be separated in space- 
hke and time-hke sectors when the spectral functions aquire a finite width. The space-hke part of 
(4) defines a potential energy density Vp since the field quanta involved are virtucll and correspond 
to partons exchanged in interaction diagrams. The time-like part of (4) corresponds to effective field 
quanta which can be propagated within the light-cone. Related separations can be made for virtuell and 
time-like parton densities [12, 13]. Without repeating the details we mention that mean-field potentials 
for partons can be defined by the derivative of the potential energy density Vp with respect to the 
time-like parton densities and effective interactions by second derivatives of Vp (cf. Section 3 in Ref. 
[13]). 

Based on the DQPM wc have developed an off-shell transport approach [14] denoted as PHSD where 
the degrees-of-freedom are dynamical quarks, antiquarks and gluons {q, q, g) with rather large masses 
and broad spectral functions in line with (1) - (3) as well as the conventional hadrons (described in the 
standard HSD approach [15]). On the partonic side the following clastic and inelastic interactions are 
included qq ^ qq, qq qq, gg ^ gg, gg ^ g. qq ^ g exploiting 'detailed-balance' with interaction 
rates from the DQPM. The hadronisation, i.e. transition from partonic to hadronic degrees of freedom, 
is described by local covariant transition rates e.g. for q + q fusion to a meson m of four-momentum 
p — {(jj, p) at space time point x — {t, x): 

X^q Pq{Pq) ^q Pq{Pq) \Vqq\^ Wm{Xq - Xq,Pg - Pg) X Ng{Xg,Pg) Ng{Xg,Pg) 5(flaV0r, Color). (7) 

In (7) we have introduced the shorthand notation Trj = J2j I d^Xjd'^pj/ {2tt)'^ where J2j denotes a 
summation over discrete quantum numbers (spin, flavor, color); Nj{x,p) is the phase-space density of 
parton J at space-time position x and four-momentum p. In Eq. (7) ^(flavor, color) stands symbolically 
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for the conservation of flavor quantum numbers as well as color neutrality of the formed hadron m which 
can be viewed as a color-dipole or 'pre-hadron'. Furthermore, Vqq{pp) is the effective quark- ant iquark 
interaction from the DQPM defined by Eq. (31) of Ref. [13] as a function of the local parton {q + q + g) 
density pp (or energy density). Furthermore, Wm{xiP) is the phase-space distribution of the formed 
'pre-hadron'. It is taken as a Gaussian in coordinate and momentum space with width v^< > = 0.66 
fm (in the rest frame) which corresponds to an average rms radius of mesons. The width in momentum 
space is fixed by the uncertainty principle, i.e. AxAp = 1 (in natural units). Related transition rates 
(to Eq. (7)) are defined for the fusion of three off-shell quarks (gi + ^2 + ^3^-8) to color neutral 
baryonic {B or B) resonances of finite width (or strings) fulfilling energy and momentum conservation 
as well as fiavor current conservation [14]. 

3 Hadronization of an expanding partonic fireball 

We now turn to actual results from PHSD for the model case of an expanding partonic fireball at 
initial temperature T — 1.7 Tc (Tc= 0.185 GeV) with quasiparticle properties and four-momentum 
distributions determined by the DQPM at temperature T = 1.7 Tg. The initial distribution for quarks, 
antiquarks and gluons in coordinate space is taken as a Gaussian ellipsoid with a spatial eccentricity 



and < >=< > 'm order to allow for the buildup of elliptic flow (as in semi-central nucleus- nucleus 
collisions at relativistic energies) . In order to match the initial off-equilibrium strange quark content in 
relativistic pp collisions the number of s (and s quarks) is assumed to be suppressed by a factor of 3 
relative to the abundance of u and d quarks and antiquarks. In this way we will be able to investigate 
additionally the question of strangeness equilibration. For more details of the numerical treatment we 
refer the reader to Ref. [14]. 

In Fig. 1 (l.h.s.) we show the energy balance for the expanding system at initial temperature 
T — l.lTc and eccentricity e =0. The total energy Etot (upper line) - which at i = is given by (4) 
integrated over space - is conserved within 3% throughout the partonic expansion and hadronization 
phase such that for t > 8 fm/c it is given essentially by the energy contribution from mesons and 
baryons (+antibaryons). The initial energy splits into the partonic interaction energy Vp (cf. Eq. (19) 
in Ref. [13]) and the energy of the time- like (propagating) partons 



with fractions determined by the DQPM [13]. In Eq. (9) the summation over i runs over all testparticles 
in an individual run. The hadronization mainly proceeds during the time interval lfm/c<t<7fm/c 
(cf. r.h.s of Fig. 1 where the time evolution of the q,q,g, meson and baryon (-|-antibaryon) number 
is displayed). As one observes from Fig. 1 on average the number of hadrons from the resonance or 
'string' decays is larger than the initial number of fusing partons. 

In order to shed some light on the hadronization process in PHSD wc display in Fig. 2 the invariant 
mass distribution of qq pairs (solid line) as well as qqq (and qqq) triples (dashed line) that lead to the 
formation of final hadronic states. In fact, the distribution for the formation of baryon (antibaryon) 
states starts above the nucleon mass and extends to high invariant mass covering the nucleon resonance 
mass region as well as the high mass continuum (which is treated by the decay of strings within the 
JETSET model [18]). On the 'pre-mesonic side the invariant-mass distribution starts roughly above the 
two-pion mass and extends up to continuum states of high invariant mass (described again in terms of 
string excitations). The low mass sector is dominated by p, ai, uj or K*, K* transitions etc. The excited 
'pre-hadronic' states decay to two or more 'pseudoscalar octet' mesons such that the number of final 
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Figure 1: (Color online) L.h.s.: Time evolution of the total energy Etot (upper solid line), the partonic 
contributions from the interaction energy Vp and the energy of time-like partons Tp in comparison to the 
energy contribution from formed mesons Em and baryons (+ antibaryons) -Eb+b- R-h.s: Time evolution 
in the parton, meson and baryon number for an expanding partonic fireball at initial temperature 
T = 1.7 Tc with initial eccentricity e = 0. 




s'" [GeV] 



Figure 2: (Color online) The invariant mass distribution for fusing qq pairs (solid line) as well as qqq 
(and qqq) triples (dashed line) that lead to the formation of final hadronic states for an expanding 
partonic fireball at initial temperature T = 1.7 with initial eccentricity e = 0. 

hadrons is larger than the initial number of fusing partons. Accordingly, the hadronization process in 
PHSD leads to an increase of the total entropy and not to a decrease as in case of coalescence models 
[6, 7]. This is a direct consequence of the finite (and rather large) dynamical quark and antiquark 
masses as well as mean-field potentials which - by energy conservation - lead to 'pre-hadron' masses 
well above those for the pseudo-scalar meson octett or the baryon octett, respectively. This solves 
the entropy problem in hadronization in a natural way and is in accordance with the second law of 
thermodynamics! 

The parton dynamics itself is governed by their propagation in the time-dependent mean-field Up{ppY 
which is adopted in the parametrized form (as a function of the parton density pp) given by Eq. (29) in 
Ref. [13]. Since the mean- field Up is repulsive the partons are accelerated during the expansion phase 
on expense of the potential energy density Vp which is given by the integral of Up over pp (cf. Sec. 3 
in Ref. [13]). The interaction rates of the partons are determined by effective cross sections which for 
gg scattering have been determined in Ref. [3] as a function of T /Tc. The latter are re-parametrized in 
the actual calculation as a function of the parton density using the available dependence of Pp{T) on 

^Analoguous mean- field phenomena are well known from studies on K^, p-meson and nucleon flow at SIS energies 
[19, 20, 21]. 
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the temperature T from the DQPM. The channels qq ^ g are described by a relativistic Breit-Wigner 
cross section which is determined by the actual masses of the fermions, the invariant energy ^/s and 
the resonance parameters of the gluon (from the DQPM). In this case a further constraint on flavor 
neutrality and open color is employed. The gluon decay to a uu, dd or ss pair is fixed by detailed 
balance. Further channels arc gg ^ g which are given by Breit-Wigner cross sections (with the gluon 
resonance parameters) and detailed balance, respectively. 

It is also interesting to have a look at the final hadronic particle ratios K~^/7t~^, p/tt"*", A/K^ etc. 
(after hadronic decays) which are shown in Table 1. The latter ratios are compared to the grand- 
canonical statistical hadronization model (SM) [22, 23] at baryon chemical potential /jLb — 0. For 
jiB = ^ the particle ratios depend only on temperature T and one may fix a freeze-out temperature, 
e.g., by the proton to tt^ ratio. A respective comparison is given also in Table 1 for T = 160 MeV and 
170 MeV for the SM which demonstrates that the results from PHSD are close to those from the SM 
for T Ri 170 MeV. This also holds roughly for the A/K+ ratio. On the other hand the K+/'k+ ratio 
only smoothly depends on the temperature T and measures the amount of strangeness equilibration. 
Recall that we initialized the system with a relative strangeness suppression factor of 1/3. The deviation 
from the SM ratio by about 13% indicates that strangeness equilibration is not fully achieved in the 
calculations. This is expected since the partons in the surface of the fireball hadronize before chemical 
equilibration may occur. 







A/K+ 


K+/n+ 


PHSD 


0.086 


0.28 


0.157 


SM T = 160 MeV 


0.073 


0.22 


0.179 


SM T = 170 MeV 


0.086 


0.26 


0.180 



Table 1: Comparison of particle ratios from PHSD with the statistical model (SM) for T= 160 MeV 
and 170 MeV. 

The agreement between the PHSD and SM results for the baryon to meson ratio in the strangeness 
S=0 and S=l sector may be explained as follows: Since the final hadron formation dominantly proceeds 
via resonance and string formation and decay - which is approximately a micro-canonical statistical 
process - the average over many hadronization events with different energy/mass and particle number 
(in the initial and final state) leads to a grand- canonical ensemble. The latter (for //^ = 0) is only 
characterized by the average energy or an associated Lagrange parameter (3 = 1/T. 

Of additional interest arc the collective properties of the partonic system during the early time 
evolution. In order to demonstrate the buildup of elliptic flow we show in Fig. 3 (l.h.s.) the time 
evolution of 

V2^{ipl-Pl)/{PI+Pl)) (10) 

for partons (solid line), mesons (long dashed line) and baryons (dot-dashed line) for an initial eccentricity 
e = 0.33. As seen from Fig. 3 the partonic flow develops within 2 fm/c and the hadrons essentially 
pick up the collective flow from the accelerated partons. The hadron V2 is smaller than the maximal 
parton V2 since by parton fusion the average V2 reduces and a fraction of hadrons is formed early at 
the surface of the flreball without a strong acceleration before hadronization. It is important to point 
out that in PHSD the eUiptic flow of partons predominantly stems from the gradients of the repulsive 
parton mean-flelds (from the DQPM) at high parton (energy) density. To demonstrate this statement 
we show in Fig. 3 (r.h.s.) the result of a simulation without elastic partonic rescattering processes by 
the short dashed line. 

We note in passing that the final hadron V2 increases linearly with the initial eccentricity e and 
indicates that the ratio V2/e is practically constant (~ 0.2) [14] as in ideal hydrodynamics. Accordingly 
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Figure 3: (Color online) L.h.s.: Time evolution of the elliptic flow V2 for partons and hadrons for the 
initial spatial eccentricity e = 0.33 for an expanding partonic fireball at initial temperature T = 1.7 T^.. 
R.h.s.: The elliptic flow V2 - scaled by the number of constituent quarks Uq - versus the transverse 
kinetic energy divided by Uq for mesons (dashed line) and baryons (solid line). 



the parton dynamics in PHSD are close to ideal hydrodynamics. This result is expected since the ratio 
of the shear viscosity 77 to the entropy density s in the DQPM is on the level of r^/s ~ 0.2 [3] and thus 
rather close to the lower bound of 77/5 = l/(47r). 

A further test of the PHSD hadronization approach is provided by the 'constituent quark number 
scaling' of the elliptic flow f 2 which has been observed experimentally in central Au+Au collisions at 
RHIC [4, 24]. In this respect we plot v^jnq versus the transverse kinetic energy per constituent parton, 

vht — m 

Tkin = , (11) 

riq 

with mj- and m denoting the transverse mass and actual mass, respectively. For mesons we have Uq = 2 
and for baryons/antibaryons Uq = 3. The results for the scaled elliptic flow are shown in Fig. 3 (r.h.s.) 
for mesons and baryons and suggest an approximate scaling. We note that the scaled hadron elliptic flow 
V2/nq does not reflect the parton V2 at hadronization and is significantly smaller. Due to the limited 
statistics especially in the baryonic sector with increasing pt this issue will have to be re-addressed 
with high statistics in the actual heavy-ion case where the very early parton pt distribution also shows 
'power-law' tails. 



4 Summary 

In summary, the expansion dynamics of an anisotropic partonic fireball has been studied within the 
PHSD approach which includes dynamical local transition rates from partons to hadrons (7) (and 
vice versa). It shows collective features as expected from ideal hydrodynamics in case of strongly 
interacting systems. The hadronization process conserves four-momentum and all flavor currents and 
slightly increases the total entropy since the 'fusion' of rather massive partons dominantly leads to the 
formation of color neutral strings or resonances that decay microcanonically to lower mass hadrons. 
This solves the entropy problem associated with the simple coalescence model! 

We find that the hadron abundancies and baryon to meson ratios are compatible with those from the 
statistical hadronization model [22, 23] - which describes well particle ratios from AGS to RHIC energies 
- at a freeze-out temperature of about 170 MeV. Our calculations show that the hadron elliptic flow 
is essentially produced in the early partonic stage where also the strong repulsive parton mean-fields 
contribute to a large extent. 
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